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ABSTRACT 



43 

Of 

We report the first detection of chloronium, HiQ + , in the interstellar medium, using the HIFI instrument aboard the Herschel Space Observatory. 
O ' The 2 l2 - loi lines of ortho-H^Cl* and ortho-H^CP are detected in absorption towards NGC 63341, and the In - Ooo transition of para-H 35 Cl + 
^3 , is detected in absorption towards NGC 63341 and Sgr B2(S). The H2C1 + column densities are compared to those of the chemically-related species 
. HC1. The derived HC1/H 2 C1 + column density ratios, ~ 1—10, are within the range predicted by models of diffuse and dense Photon Dominated 
, ^ , ■ Regions (PDRs). However, the observed HtCP column densities, in excess of 10 13 cm -2 , are significantly higher than the model predictions. Our 
observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints 
I , for astrochemical models. 

^ [ Key words. Astrochemistry — ISM: abundances — ISM: molecules — Line: identification — Molecular processes — Submillimetre: ISM 

^-H 1 ■ Introduction The product of this reaction is the reactive HC1 + ion, which un- 

r • ' dergoes further reaction with H2 to form H2CL: 

L The halogen elements, fluorine and chlorine, form hydrides 

that are very strongly bound: hydrogen fluoride is the only di- HC1 + + H2 — > H2CL + H. (2) 

atomic hydride, and HC1 + the only diatomic hydride cation, 

with a dissociation energy exceeding that of molecular hydro- The H 2C1 + molecule does not react with H 2 , and is destroyed 
gen. Drawing upon earlier work by Jura (1974), Dalgarno et b Y dissociative recombination and proton transfer to CO, both of 
al. (1974), van Dishoeck & Black (1986), Blake et al. (1986), which are sources of hydrogen chloride, HC1 
Schilke et al. (1995), Federman et al. (1995), and Amin (1996), Prior t0 the launch of Herschel, the H CI and H CI iso- 
b ■ Neufeld & Wolfire (2009; hereafter NW09) have recently car- topologues were the only chlorine-containing molecules to have 
ried out a theoretical study of the chemistry of chlorine-bearing been detected in the ISM (e.g., Blake et al. 1985; Zmuidzinas 
molecules, in both diffuse and dense molecular clouds. In diffuse et al - 1995 ; Schilke et al. 1995; Salez et al. 1996; see also re- 
interstellar gas clouds, the dominant ionization state of every el- cent HIFI observations of Cernicharo et al. 2010)U However, 
ement is determined by its ionization potential. Chlorine, with an predictions for the chemistry of Cl-beanng interstellar molecules 
ionization potential slightly lower than that of hydrogen, is pre- (NW09) have identified chloronium, H 2 C1 + , as a relatively abun- 
dominantly singly-ionized. The Cl + ion can react exothermically dant s P ecies that 18 potentially detectable. H 2 C1 + is predicted to 
with H 2 , the dominant molecular constituent of the interstellar be most abundant in those environments where the ultraviolet ra- 
medium (ISM)' diation is strong: in diffuse clouds, or near the surfaces of dense 

clouds that are illuminated by nearby O and B stars. In such 

„ 1+ j j „„+ j j , 1N environments, the photoionization of atomic chlorine leads to 

CI + H2 — > HC1 + rl. (1) 



o 



X 



1 The metal halides NaCl, KC1, and A1C1 have been detected in the 

* Herschel is an ESA space observatory with science instruments circumstellar envelope of the evolved star IRC+ 10216, with abundances 

provided by European-led Principal Investigator consortia and with im- that reflect the thermochemical equilibrium established within the stel- 

portant participation from NASA. lar photosphere (Cernicharo & Guelin 1987). 
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a large abundance of Cl + ions that can form HC1 + and H2C1 + 
through reactions (1) and (2). A secondary abundance peak oc- 
curs in dense, shielded regions; here HQ becomes a significant 
reservoir of gas-phase chlorine, and can produce H2C1 + through 
reaction with H^ : 

HC1 + H+ -> H 2 C1 + + H 2 . (3) 

However, the chlorine depletion is typically large within such re- 
gions (Schilke et al. 1995) and thus the overall H2C1 + abundance 
is rather small. 

In diffuse molecular clouds of density «h = 10 2 ' 5 cm" 3 , H 2 
column density > 10 20 cm" 2 , and^uv in the range 1-10 (where 
Xuv is the UV radiation field normalized with respect to the 
mean interstellar value, Draine 1978), the NW09 model pre- 
dicts H 2 CF column densities ~3xl0 10 yuv cm" 2 . In dense PDRs 
(«h = 10 4 cm" 3 ) illuminated by strong radiation fields (xuv > 
10 3 ), the predicted H 2 C1 + column densities are ~10 12 cm" 2 . 

In this Letter, we report the first detection of chloronium to- 
wards NGC 63341 and Sgr B2(S), obtained using the HIFI in- 
strument (de Graauw et al. 2010) aboard the Herschel Space 
Observatory (Pilbratt et al. 2010). NGC 6334 is a luminous and 
relatively nearby (1.7 kpc) molecular cloud/Hn region complex 
containing several concentrations of massive stars at various 
stages of evolution. The far-infrared source "I", located at the 
northeastern end of the complex, is associated with a NIR clus- 
ter of bolometric luminosity of 2.6 x 10 5 L (Sandell 2000), with 
four embedded compact millimeter continuum sources (Hunter 
et al. 2006). Sgr B2(S) is a strong submillimeter continuum 
source with a much less complex hot core emission spectrum, as 
compared to its better known neighbor Sgr B2(M). This makes 
it a prime candidate for absorption studies, probing the entire 
sight-line between the Sun and the Galactic center, with clouds 
in the Orion, Sagittarius, and Scutum spiral arms easily identi- 
fied at separate velocities (e.g., Greaves & Nyman 1996). 

2. Observations 

HIFI observations presented here were carried out between 2010 
March 1 and March 23, using the dual beam switch (DBS) ob- 
serving mode, as part of guaranteed and open time key pro- 
grams CHESS: Chemical HErschel Spectral Surveys, HEXOS: 
Herschel/HIFI Observations of EXtraOrdinary Sources: The 
Orion and Sagittarius B2 starforming regions, and HOP: 
Herschel Oxygen Program. The source coordinates are: a y 2 ooo = 
17 h 20 m 53.32 s and 6 j20 qo = -35°46'58.5" for NGC 63341, and 
a/2000 = 17 h 47 m 20.3 s and5 y20 oo = -28°23'43.0" forSgrB2(S). 
The DBS reference beams lie approximately 3' east and west 
(i.e. perpendicular to the roughly north-south elongation of the 
two sources). Because the DBS mode alternates between two 
reference positions, separated by 6' on the sky, we used the 
Level 1 data to compute a difference spectrum between the 
two reference positions to check for possible contamination in 
the reference beams; we see no evidence for emission or ab- 
sorption in such a difference spectrum. We used the HIFI wide 
band spectrometer (WBS) providing a spectral resolution of 
1.1 MHz (-0.4 kms" 1 at 780 GHz) over a 4 GHz IF band- 
width. The spectra presented here are averages of the H and V 
polarizations, with equal weighting, reduced using HIPE (Ott 
2010) with pipeline version 2.6. The resulting Level 2 DSB 
spectra were exported to the FITS format for a subsequent 
data reduction and analysis using the IRAM GILDAS package 
( http://iram.fr/IRAMFR/GILDAS). 

The band 2b, lb and la spectral scans of NGC 63341 consist 
of double sideband spectra (DSB) with a redundancy of 8, which 



gives observations of a lower or upper sideband frequency with 8 
different settings of the local oscillator (LO). The Sgr B2(S) data 
consist of 8 LO settings with a high redundancy of 12, centered 
near the frequency of the 487.2 GHz line of 2 . The observa- 
tions were fine-tuned so that 4 of the 8 LO settings cover the 
frequency of the p-H 2 Cl + line. These observing modes allow for 
the deconvolution and isolation of a single sideband spectrum 
(Comito & Schilke 2002). We applied the standard deconvolu- 
tion routine within CLASS. All NGC 63341 data presented here 
are deconvolved single sideband spectra, including the contin- 
uum. The HC1 data in Sgr B2(S) were obtained using the DBS 
single point observing mode with 3 shifted LO settings that were 
averaged to produce the final spectrum. The HIFI beam size at 
485 GHz and 780 GHz is 44" and 30", respectively, with main 
beam efficiency of ~0.68. 

3. Spectroscopy of H 2 CI + 

The H 2 C1 + ion is a closed-shell molecule, isoelectronic with 
H 2 S. Like H 2 S and water, H 2 C1 + is a highly asymmetric top, ex- 
hibiting a b-type rotational spectrum. Its fairly large dipole mo- 
ment, calculated ab initio to be 1 .89 D (Miiller et. al 2005), about 
70% larger than that of HC1 (1.109 D; de Leluw & Dymanus 
1973), results in strong lines in the THz range. Araki et al. (2001) 
measured rotational spectra of H 35 CF, H 2 37 C1 + , and HDC1 + be- 
low 500 GHz The accurate spectroscopic constants derived from 
these measurements, including electric quadrupole coupling pa- 
rameters, yield a central bond angle in H 2 C1 + of ~94.2° (sim- 
ilar to that of H 2 S, 92.2°; Burrus & Gordy 1953), and permit 
the prediction of the ground state ortho transitions of H 2 C1 + and 
H 2 37 C1 + near 780 GHz to well within 1 MHz (see Table 1 in the 
Online Material). 

4. Results 

4.1. NGC 63341 

Strong absorption at the frequency of the 2j 2 - loi transition of 
o-H 35 Cl + at 781.6 GHz in NGC 63341 (Fig. la) has provided 
the initial identification. Fitting the o-H^CF 2i 2 - Iqi hyper- 
fine structure (HFS) gives a line velocity of -1.7 kms" 1 for the 
strongest hyperfine component and a line width of 11.6 kms" 1 . 
The OH absorption profiles (Brooks & Whiteoak 2001) re- 
veal two molecular clouds located along the line of sight to 
NGC 6334, one with velocities extending from -15 to 2 kms" 1 , 
and one with a well-defined velocity near 6 kms" 1 . The hot 
core emission lines peak at about -6.5 kms" 1 (e.g., C ls O 7- 
6 in the same band; HC1, Sect. 4.1). Water and CH spectra 
towards NGC 63341 show multiple velocity components, in- 
cluding absorption near kms" 1 , close to the H| Cl + veloc- 
ity (Emprechtinger et al. 2010; van der Wiel et al. 2010). The 
chloronium line velocity in NGC 63341 is in good agreement 
with the H 2 + absorption velocity, also a tracer of diffuse gas, 
when H 2 + frequencies of Miirtz et al. (1998) are used (see 
Schilke et al. 2010 for a discussion of the H 2 + line frequen- 
cies). The large H 2 C1 + line width may be due to blending of 
multiple absorption components. However, the H 2 + line width 
is also quite large, about 8 km s" 1 . 

The corresponding line of o-H?, 7 Cl + is also detected 
(Fig. lb). However, the spectrum is contaminated by interfering 
emission of dimethyl ether (light-blue line in Fig. lb), one of 
the most abundant "weeds" in NGC 63341. The contamination 
is subtracted by using an LTE model that fits profiles of nearby 
dimethyl ether lines with similar upper level energies (Endres et 
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Fig. 1. Spectra of chlorine species in NGC 63341: (a) o-H 35 CF 
2 12 - lot, (b) o-H 37 CF 2i 2 - lot, (c) p-H 35 CF 1„ - 00 , (d) 
H 35 C1 1-0, and (e) H 37 C1 1-0. The velocity scale corresponds 
to the strongest HFS components. Green lines show HFS fits and 
positions of the HFS components. The o-H?, 7 CF and p-H?, 5 CF 
lines are blended with dimethyl ether emission (light-blue lines 
in panels b and c). 



al. 2009). The resulting o-H^CF spectrum is shown as a black 
line in Fig. lb. The In - Ooo line of p-H| 5 Cl + (Fig. lc) is also 
blended with dimethyl ether emission, similarly subtracted. 



FbCF spectra in NGC 63341 can be compared to those of 
the chemically related species HC1 (Fig. Id and le). The lines 
of H 35 C1 and H 37 C1 are detected in emission at the hot core ve- 
locity (~-6.3 kms~'), with narrow line widths of 4.1 and 3.3 
kms , respectively. The HC1 HFS is resolved spectrally using 
HIFI WBS, allowing for the determination of the line opacity. 

We have modelled the H2CF spectra assuming the same ex- 
citation temperature of 5 K for all hyperfine components A low 
value of the excitation temperature is justified given the high 
spontaneous emisison rates and critical densities of the transi- 
tions considered here and it provides a lower limit for the molec- 
ular column densities derived from absorption measurements. 
An HFS fit to the 2\2 - lot transition of o-H^CF (green line 
in Fig. la) gives an o-H^CF column density of 1.3 x 10 13 cm -2 , 
under the assumption that the absorption completely covers the 
continuum and is not concentrated in small clumps. A fit to 
the In - Ooo spectrum of p-H^CF, with all parameters other 
than the column density fixed, gives a p-H^CF column den- 
sity of 4.0 x 10 12 cirT 2 . The total H 3 , Cl + column density is thus 
1.7 x 10 13 cirT 2 and the ortho-to-para ratio is 3.2, consistent 
with the statistical weight ratio. For an excitation temperature 
of 2.7 K, the ortho and para H2CF column densities are approx- 
imately 10% and 20% lower, respectively. The H^CF spectra 
are all optically thin (line center optical depth of ~0.2 for the 
ortho line). We derive an H^CF/H^CF ratio of 3, close to the 
terrestrial ratio of 3. 1 . 

We have modelled the H 35 C1 and H 37 C1 emission spectra as- 
suming a source size of 10" (approximate size of the cluster of 
compact continuum sources seen in the SMA image of Hunter et 
al. 2006). Under this assumption, a least squares fit to the H 35 C1 
spectrum gives an excitation temperature of 3 1 K and a column 
density of 4.0 x 10 14 cirT 2 . For H 37 C1, we derive an excitation 
temperature of 2 1 K and a column density of 1 .5 x 1 14 cirT 2 ; the 
resulting H 35 C1/H 37 C1 ratio is 2.7. However, HC1 column densi- 
ties and the isotopic ratio depend strongly on the assumed source 
size (for a source size of 5" the derived isotopic ratio is 4. 1). The 
350 jum continuum flux density toward NGC 63341 is 1430 Jy 
in a 9" beam (CSO/SHARC II; Dowell et al., private comm.) 
Assuming a dust temperature of 100 K (Sandell 2000) and a 
grain emissivity K350 = 0.1 cm~ 2 g~', we derive an H2 column 
density of 1.2 x 10 24 cirT 2 , which implies an H 35 C1 abundance 
of -1.7 x 10" 10 with respect to H nuclei. The lines of both HC1 
isotopologues are optically thick, with line center optical depths 
of ~2.2 and 1.6 for the strongest hyperfine component of H 35 C1 
and H 37 C1, respectively. 

4.2. Sgr B2(S) 

The p-H^CF spectrum towards Sgr B2(S) (Fig. 2, upper panel) 
shows strong absorption near the systemic velocity of the Sgr B2 
envelope (~62 km s~') and two additional deep absorption com- 
ponents between and 20 kms 1 . In addition, shallow ab- 
sorption is seen over a broad range of velocities down to 
-100 kms 1 , in agreement with the Hi absorption spectrum to- 
wards the nearby source Sgr B2(M) (magenta line in Fig. 2). 
Both H 35 C1 and H 37 C1 (Fig. 2, lower panel) show deep absorp- 
tion at the envelope velocity and a shallow absorption between 



2 We made use of the myXCLASS program 
( http://www.astro.uni-koeln.de/projects/schilke/XCLASS ), 
which accesses ~the CD MS (Miiller at al. 2001, Miiller et 
al. 2005; |http://www.cdms.de) and JPL (Pickett et al. 1998; 
http://spec.jpl.nasa.gov I molecular databases. 
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Fig. 2. Spectra ofp-H 35 Ci + lii-0 o (upper panel), andH 35 Cl and 
H 37 C1 1 - (lower panel) towards Sgr B2(S). The H 37 C1 spec- 
trum has been shifted up by 0.5 K. The magenta line in the up- 
per panel shows the H i absorption spectrum towards Sgr B2(M) 
(from Garwood & Dickey 1989). 



and 20 kms -1 . Similarly to NGC 63341, we see velocity offsets 
of order a few kms 1 between HQ and H2C1 + components. 

Assuming a 5 K excitation temperature (the same as for 
NGC 63341) and an ortho/para ratio of 3, we derive H2C1 + col- 
umn densities of 3.4 x 10 13 and 2.2 x 10 13 cirT 2 for the and 
62 kms 1 components, with corresponding H 35 C1 column den- 
sities of 4 x 10 13 and 2 x 10 14 cm 2 . The H 35 C1/H 37 C1 ratio is 
~3.3 in both components. We estimate the uncertainties in our 
molecular column density estimates to be of order a factor of 2. 

To derive the hydrogen column density in the foreground gas 
towards Sgr B2(S), we use the method employed in Lis et al. 
(2001) to analyze the Oi absorption towards Sgr B2(M), based 
on H i and 13 CO absorption data. We assume that the foreground 
absorption is extended and column densities are the same to- 
wards Sgr B2(M) and (S). We derive a total hydrogen nuclei 
column density of ~2 x 10 22 cm -2 in the atomic and molecular 
components in the velocity range -10 to 20 kms -1 (with a fac- 
tor of 2 uncertainty). The corresponding chlorine content, in the 
form of H 2 C1 + and HC1, is 7 x 10 13 cm 2 , implying a Cl/H ratio 
of ~4x 10 -9 . This can be compared to the values measured in the 
UV in diffuse clouds (e.g., Sonnentrucker et al. 2006), which are 
in the range 3x 10~ 8 -4x 10~ 7 . Therefore the high H 2 C1 + column 
densities we derive here are consistent with the overall chlorine 
budget, leaving plenty of room for atomic CI and depletion on 
dust grains. 



5. Discussion 

Our estimates of the H2CF column densities towards 
NGC 63341 and Sgr B2(S), in excess of 10 13 cirT 2 , are signif- 



icantly higher than those expected for a single dense or diffuse 
PDR viewed at normal incidence. This might point to some defi- 
ciency in the models. Alternatively, a significant enhancement in 
the absorbing column density could result if the normal to the ir- 
radiated surface were inclined relative to the sight-line, or indeed 
if multiple PDRs were present along the sight-line, particularly if 
the radiation field is enhanced, as may be likely for the multiple 
absorption components seen towards Sgr B2. Similar discrepan- 
cies between models and observations are seen for other reactive 
ions in massive starforming regions (e.g., CO + toward AFGL 
259 1 ; Bruderer et al. 2009). 

We derive an HC1/H 2 C1 + ratio of -10 in NGC 63341 and 
the Sgr B2 envelope (assuming that in the case of NGC 63341 
the H2CF column density on the back side is the same as 
that derived in front of the continuum source from our absorp- 
tion measurements). This is well within the range predicted 
for dense PDRs (up to ~100 for densities above 10 6 cm -3 ). 
The HCI/H2CL ratio derived in the foreground gas towards 
SgrB2(S) at velocities 0-20 kms -1 , ~1, is also consistent with 
predictions of diffuse cloud models. 

While a detailed analysis of chlorine chemistry in these and 
other sources that have been or will be observed using HIFI will 
be presented in a forthcoming paper, this work clearly demon- 
strates the outstanding spectroscopic capabilities of HIFI in the 
search for new interstellar molecules, particularly hydrides, and 
in providing robust constraints for astrochemical models of the 
interstellar medium. 
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Table 1. Frequencies of the H2CF transitions observed. 





Transition 


Frequency 


Error 




Ei 








( MHz) 


(MHz) 


(s ') 


/ -i\ 
(cm ') 








H^ 5 C1 + 








Ill 


- Ooo 3/2 - 


3/2 


485413.427 


0.029 


0.00159 







5/2- 


3/2 


485417.670 


0.015 


0.00159 







1/2- 


3/2 


485420.796 


0.057 


0.00159 





2l2 


- loi 5/2- 


5/2 


781609.303 


0.063 


0.00179 


14.1 




3/2- 


1/2 


781611.062 


0.062 


0.00248 


14.1 




5/2- 


3/2 


70 1 tin 71 1 

/5I0ZZ. /zl 


0.063 


A AA A 1 H 

0.00417 


14.1 




7/2- 


5/2 


781626.794 


0.060 


0.00596 


14.1 




1/2- 


1/2 


781628.554 


a i 

0.061 


A AA/1A<; 

0.00496 


i/ii 
14.1 




3/2- 


3/2 


781635.214 


a a<;o 

0.062 


A AA1 1 O 

0.00318 


i/ii 
14.1 








HfCF 








2l2 


- loi 5/2- 


5/2 


780037.315 


0.069 


0.00178 


14.1 




3/2- 


1/2 


780038.760 


0.066 


0.00247 


14.1 




5/2- 


3/2 


780047.903 


0.066 


0.00414 


14.1 




7/2- 


5/2 


780051.197 


0.062 


0.00592 


14.1 




1/2- 


1/2 


780052.642 


0.065 


0.00493 


14.1 




3/2- 


3/2 


780057.820 


0.068 


0.00316 


14.1 



Note: Frequencies and spontaneous emission coefficients have been 
calculated from the constants derived by Araki et al. (2001), see also 



CDMS. 
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